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It is demonstrated analyti(;ally that the capillary depression C of a nearly planar meniscus is linearly 
related to its meniscus lieijiht h by the relationship C = /i/[l —/()(r/a)]. Here, a^=ylpg, p and 7 are 
the density and surface tension of the licpiid respectively,^, the acceleration due to jiravity, r is the radius 
of the capillary tube at th<' licpiid-solid-^as interface, and L){z) is the modified Bessel function of the 
first kind and zero order. 
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1 . Introduction 

Corrections due to capillarity must be known or 
eliminated in accurate manometric procedures. The 
problem of capillary depression (or rise) for arbitrarily 
shaped gas-liquid interfaces has been studied ex- 
tensively, both theoretically and experimentally. Little 
attention, however, has been given specifically to 
nearly planar menisci occurring in large bore precision 
mercury manometers and which have been produced 
by Moser and Poltz [1] ^ in the compression capillary 
of a McLeod gage to reduce errors due to capillarity. 

It is well known [2, 3, 4] that a finite pressure differ- 
ence exists across a gas-liquid interface whenever it 
is not planar. Laplace [2] was the first to show theo- 
retically that such local pressure differences were 
solely determined from the mean curvature of the sur- 
face and the surface tension of the liquid. This is a 
statement of Laplace's well-known nonlinear differ- 
ential equation for liquid manometric work. Numerical 
solutions [3, 4, 5, 6] of Laplace's equation exist and 
tables showing the dependence of capillary depres- 
sions on meniscus heights and tube bores have been 
constructed. Values appearing in these tables are for 
arbitrary menisci, and interpolation is frequently 
required when using them. 

In this work only nearly planar menisci are con- 
sidered, and their associated capillary depressions are 
derived analytically. 

2. Theory 

The method of the calculus of variations [7] has been 
used by Goodrich [8j in his study of the static meniscus. 
This treatment of the static meniscus is similar to that 



of Goodrich in that the calculus of variations is apphed. 
Figures la and lb represent the cross section of 
portions of capillary systems which are capable of 
producing nearly planar menisci which are of capillary 
dimensions [1]. The coordinate systems and variables 
used in this discussion are shown in these figures. In 
these figures LMN and L'M'N' represent portions of 
the vessel. MO'M' represents the hquid-gas interface. 
OO'O" is taken to be the axis of symmetry of the 
meniscus. The ;ty-plane is perpendicular to OO'O" with 
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' Figures in brackets indicate the literature references at the end of this paper. 



Figure L Coordinate system and displacements for a nearly 
planar meniscus. 
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zix, y) = \ x^ + j2 and ^= tan"^y/jc (not shown in these 
figures). This distance from 00' to the gas-liquid- 
sohd interface is r. AA' represents the position the 
meniscus would take in the absence of capillarity. 
All displacements i//(jc, y) are taken with respect to 
AA', 

Let Wy and wy be energies associated with the actions 
of gravity and surface tension, respectively. Consider 
the meniscus generated from the rotation of the seg- 
ment O'M' about 00'. For convenience of discussion, 
consider the situation shown in figure lb. The potential 
energy of a column of liquid of unit cross-sectional 
area and height i//(jc, y) is given by 



Wfj = i pgijj^ix, y). 



(1) 



where p is the density of the liquid, g is the accelera- 
tion due to gravity. The surface element zdOdt^, where 
dt, is an arc length measured along O'M', possesses a 
surface free energy of amount yzdOd^, where y is the 
surface tension of the liquid. This exceeds that of a 
planar element zdOdz by amount Wy = yzd6{d^ — dz) . 
Elementary geometrical considerations give 



dc = dz vTTJWP, 

where Vi//= {dijj/dx^ dipldy). Thus the excess surface 



free energy per unit area, i.e., zdSdz— 1 is given by 



term yVi//-Vx may be replaced by use of the identity 

V • (xVi//) = Vx • Vi;/+ xV^i//. Furthermore, 

V'(xVi//)=— -h — , 
dx dy 

where P = X^y-> Q = X^j- and i//.,-, ifj^ are partial deriva- 
tives of i// with respect to x and y. Then by Stokes' 
theorem 

eq (5) becomes 

dW=el I {pg\\}-yV^)xdxdy^-e \ y(t)oxds 

+ e I yxi^xdy-ifjydx). 

In terms of the normal and tangential directions at F, 
the differential elements dx, dy and the partial deriva- 
tives i//j-, ijjij become —sin Ods, cos Sds, 

- — cos ^ and:^ — sin t^, 
bn dn 

respectively, so that ^jdy — ifjydx can be written as 



U^' 



(z) dn 



-y[l- Vl-h|Vi//|2]. 

The total energy of the static meniscus is taken to be 8W= e I I {pg\\f —y\/'^\\})xdxdy 

W[x\f{x,y)]=\ {Wfj+wy)dxdy^\ (\>{s)ds. (3) C / Si1j\ , 

J J'^ J' -\-ey\Uo-^j^jxds. (6) 

The second term on the right is a hue integral term and 

has been included to take into account any energy The necessary condition that ^ is a minimum is 6^= 

associated with the hquid-gas-sofid interface along the for arbitrary variations x(^, y) [7]. Hence, we obtain 

curve F : x^-\-y^= r^. The double integral is to be taken the equations 

over the surface of the meniscus R. For |Ai//| <^ 1, we , 

have VHix.y) Mx.y)=^ (7) 



W[_^\^{x,y)]=\^ ^^{pg^\^^^y\y^\s\^)dxdy 

-f \ <i>{s)ds. 



and 



V2i//(x,y) -^(^,y)=0 



dn 



+ (^o=-0, 



(8) 



^^ where a^ = y/pg. Solutions of eq (7) give functions that 
minimize W, while eq (8) provides the condition for the 
Suppose i//*(x, y) and (/)*(5) are chosen such that evaluation of the constants of integration. It is eq (7) 



i//*(x, y) = i//(^, y)-hex{x, y) 

and (/)*(5) =(/)(5) +ey(/)()(5)x(5), where e is a very 
small number. Then the variation in the total energy, 
dW= r[i//*(x, y) ] - r[i//U, y) ] , is given by 



which is now pertinent. 

Since the gas-Hquid interface is taken to be a surface 
of revolution, i// is a function of z alone and eq (7) 
transformed to cylindrical coordinates (z, 6) becomes 



d^i\f{zla) 1 d\\f{zla) 



8W=e 



I I {pg\ljx-^y'^4^-Vx)dxdy 



dz'- 



- + 



dz 



i//U/a)=0 



(9) 



+ 0(6'^), 



(5) 



where O(e^) represents the term of order e^. The 



A solution of eq (9) is 

ilj{zla)=CJ,iz/a) (10) 

where loiz/a) is the modified Bessel function of the 
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first kind and zero order. The meniscus height h is 
given by tlie difference of ifj at z=0 and z=r, and 
eq (10) becomes 

hIo{z/a) 



\\}{zla)=-. 



(11) 



l-/()(r/a) 

where the plus and minus signs refer to ca[)illary 
depression and rise, respectively. The capillary 
depression C = i//(0) (see figure 1) is Hnearly propor- 
tional to the meniscus height, h for a specific vahic of 
the ratio ria. Differentiation of eq (10) gives 

hl\ {zja) 



V^{zla)=-. 



(12) 



a[l-/o(r/a)] 

where I\{z) is the modified Bessel function of the first 
kind and order one. The condition for validity of linear 
behavior between C and h is that | \/\\f{rla) \ < 1. 

3. Results and Discussion 

In eq (10), vertical displacements of the meniscus 
are taken with respect to a carefully selected planar 
surface such that the value of the capillary depression 
is given by C. The condition that a meniscus be nearly 
planar implies that h < r. The condition that only 
menisci of capillary dimensions are to be considered 
requires that < rja < 1, since a ~ 2 mm for mercury. 
However, larger values of rja are admissible so long 
as the condition that the meniscus be nearly planar 
is met. 

Table 1 gives capillary depressions C = i//(0), from 
eq (11) for some values of rja. Table 2 gives a com- 
parison of the capillary depressions predicted by 
eq (11) to those obtained from numerical solutions of 
Laplace's equation as performed by Akiyama, 
Hashimoto, and Nakayama [5], and (/ould and Vickers 
[6]. Table 2 shows that results of this paper are in fair 
agreement with those obtained by numerical methods. 

Figure 2 shows the dependence of capillary 
depressions on meniscus heights for some values 
of rja. The solid lines are the linear behavior pre- 

Table 1. Rdtio of the capilliiry depression to the meniscus height 
of a nearly planar mercury meniscus for some values o/r/a 



rIa 


.20 


.30 


.40 


.50 


.60 


.70 


.80 


.90 






Ratio 


100 


44.2 


24.7 


15.7 


10.9 


7.91 


6.01 


4.69 







Table 2. Comparison of capillary depressions for a meniscus height 
of mercury equal to 0.1 mm 



rIa 


AHN* 


GV** 


This work 


Difference 


|V^| 


.34 


mm //A' 

3.34 


mm H^ 


mm Hg 

3.43 
1.33 
.578 
.409 


Percent 

2.9 
0.8 
0.9 
2.0 


.3 


.543 


1.32 
.573 


.2 


815 




.1 


96 


.417 


.1 









dieted by this paper. The dotted lines are loci of 
constant Vi//. These dotted lines api)roach the axes 
asymptotically as r/a approaches zero or becomes 
very large. The points represent values obtained from 
numerical solutions of Laplace's equation. This figure 
shows that within the validity of the assumptions, 
capillary depressions are linearly related to meniscus 
heights. The figure shows that this linear relationship 
compares fairly well with results obtained from numer- 
ical solutions of Laplace's equation for | Vi// | < 0.5 
and r/a > 0.35. 
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Fl(;URE 2. Capillary depression versus meniscus height for nearly 
planar mercury menisci. 



Open and solid circles indicate values from numerical solutions of Laplace's equatic 
open circles, Gould and Vickers; solid circles, Akiyama, et al. Dotted lines indicate valu 

of I V(M. 



alues 



In summary, it has been shown that capillary 
depressions for nearly planar menisci are linearly 
related to their meniscus heights. Conditions for the 
validity of this behavior have been given. 



*Values taken from Akiyama, et al. [5]. 
**Values taken from Gould and Vickers [6]. 



The author is indebted to Stanley Ruthberg for 
suggesting this i)r()bleiTi. The author is grateful to 
H. P. Waters for some helpful discussions in connection 
with this work. 
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